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O ' ABSTRACT 

o . 

I The two-fluid effects on tlie radial outflow of relativistic electron-positron plasma 

• are considered. It is shown that for large enough Michel (1969) magnetization param- 

I eter cr 3> 1 and multiplication parameter A = n/ncj ^ 1 one-fluid MHD approxima- 

■ tion remains correct in the whole region \E\ < \B\. In the case when the longitudinal 

electric current is smaller than the Goldreich-Julian one, the acceleration of particles 

■ near the light surface \E\ = \B\ is determined. It is shown that, as in the previously 
\ considered (Beskin Gurevich & Istomin 1983) cylindrical geometry, almost all electro- 

^ ' magnetic energy is transformed into the energy of particles in the narrow boundary 

lO ' layer Avd/zj A"""^. 
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1 INTRODUCTION 



(-H [ Despite the fact that the structure of the magnetosphere of radio pulsars remains one of the fundamental astrophysical 
rS i' problems, the common view on the key theoretical question - what is the physical nature of the neutron star braking - is 
^ I absent (Michel 1991, Beskin Gurevich & Istomin 1993, Mestel 1999). Nevertheless, very extensive theoretical studies in the 
5_j ' seventies and the eighties allowed to obtain some model-independent results. One of them is the absence of magnetodipole 
[ energy loss. This result was first obtained theoretically (Henriksen & Norton 1975, Beskin et al 1983). It was shown that the 
' electric charges filling the magnetosphere screen fully the magnetodipole radiation of a neutron star for an arbitrary inclination 
^ ' angle x between the rotational and magnetic axes if there are no longitudinal currents fiowing in the magnetosphere. Later 
• this result was also confirmed by observations. The direct detections of the interaction of the pulsar wind with a companion 

star in close binaries (see e.g. Djorgovsky & Evans 1988, Kulkarni & Hester 1988) have shown that it is impossible to explain 
• the heating of the companion by a low-frequency magnetodipole wave. 

On the other hand, the detailed mechanism of particle acceleration remains unclear. Indeed, a very high magnetization 
parameter a (Michel 1969) in the pulsar magnetosphere demonstrates that within the light cylinder r < i?L ~ c/Q the main 
part of the energy is transported by the Poynting fiux. It means that the additional mechanism of particle acceleration must 
work in the vicinity of the light cylinder. It is necessary to stress that an effective particle acceleration can only take place 
for small enough longitudinal electric currents / < Igj when the plasma has no possibility to pass smoothly through the 
fast magnetosonic surface and when the light surface \E\ = \B\ is located at a finite distance. As to the case of the large 
longitudinal currents / > Igj, both analytical (Tomimatsu 1994, Begelman & Li 1994, Beskin et al 1998) and numerical 
(Bogovalov 1997) considerations demonstrate that the acceleration becomes ineffective outside the fast magnetosonic surface, 
and the particle-to-Poynting fiux ratio remains smaU: ~ cr"^/^ (Michel 1969, Okamoto 1978). 

The acceleration of an electron-positron plasma near the light surface was considered by Beskin Gurevich and Istomin 
(1983) in the simple ID cylindrical geometry for / <^ Igj- It was shown that in a narrow boundary layer Azu/uj ~ 1/A almost 
all electromagnetic energy is actually converted to the particles energy. Nevertheless, cylindrical geometry does not provide 
the complete picture of particle acceleration. In particular, it was impossible to include self-consistently the disturbance of 
a poloidal magnetic field and an electric potential, the later playing the main role in the problem of the plasma acceleration 
(for more details see e.g. Mestel & Shibata 1994). Hence, a more careful 2D consideration is necessary. 
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In Sect. 2 we formulate a complete system of 2D two-fluid MHD equations describing the electron-positron outflow 

from a magnetized body with a monopole magnetic field. The presence of an exact analytical forco-frcc solution (Michel 
1973) allows us to linearize this system which results in the existence of invariants (energy and angular momentum) along 
unperturbed monopole field lines similar to the ideal one-fluid MHD flow. In Sect. 3 it is shown that for cr 1 and A 1 
(A = n/ncj is the multiplication factor) the one-fluid MHD approximation remains true in the entire region within the light 
surface. Finally, in Sect. 4 the acceleration of particles near the light surface \E\ = \B\ is considered. It is shown that, as in 
the case of cylindrical geometry, in a narrow boundary layer Aw/tu ~ A~^ almost all the electromagnetic energy is converted 
into the energy of particles. 



2 BASIC EQUATIONS 

Let us consider a stationary axisymnictric outflow of a two-component plasma in the vicinity of an active object with a 
monopole magnetic fleld. It is necessary to stress that, of course, the monopole magnetic field is a rather crude approximation 
for a pulsar magnetosphere. Nevertheless, even for a dipolc magnetic fleld near the origin, at large distances r S> J?l in the 
wind zone the magnetic field can have a monopole— like structure. For this reason the disturbance of a monopole magnetic 
field can give us an important information concerning particle acceleration far from the neutron star. 
The structure of the flow is described by the set of Maxwell's equations and the equations of motion 
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Here E and B are the electric and magnetic fields, pe and j are the charge and current densities, and and are the 
speed and momentum of particles. In the limit of infinite paxticle energy 
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is the exact solution of Maxwell's equations. In this case 
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which coincides with the well-known Michel (1973) solution. Here 7 is the Lorentz- factor of particles, Bs is the magnetic 
fleld on the surface of the sphere r = Rb, and ps = const. As a result, the angular velocity can be rewritten in a form 
Q = 27rc|ps|/Ss- The limit 7 ^ 00 just corresponds to zero particle mass in the force-free approximation. 
It is convenient to introduce the electric fleld potential $(r, 0), so that E = — V$ and 
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and the flux function "if{r,9), so that 
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Then one can seek the first-order corrections for the case w 7^ c in the following form: 
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Such a choice corresponds to a constant particle-to-magnetic flux ratio. Here A 3> 1 is the multiplication parameter (A = 
eris/lpsl, where tIs is the concentration of particles on the surface r = Rb) which is 10^ — 10^ for radio pulsars. In what follows 
we consider for simplicity the case A — const. 

Substituting (|)-(|l|) into equations (|^)-(^, we obtain to the first-order approximation the following system of equations: 
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is the Michel's (1969) magnetization parameter, m is the electron mass, and all deflecting functions are supposed to be ^ 1. 
It is necessary to stress that for applications the magnetic field Bs is to be taken near the light cylinder Rs ~ Rl because in 
the internal region of the pulsar magnetosphere B oc r~^. As it has already been mentioned, only outside the light cylinder 
the poloidal magnetic fleld may have quasi monopole structure. As a result, 



iXmc^ 



lO^BiaA-T'p- 



(30) 



where Bo - magnetic fleld on the neutron star surface r = R. Hence, for ordinary pulsars {P ~ Is, Bo ~ lO^'^G) we have 



10* 



10^ and only for fast ones (P ~ 0.1 - O.Ols, Bo ~ lO^^G) we have cr ~ 10^ - 10^ 



Formally, this system of equations requires twelve boundary conditions. We consider for simplicity the case Q,R/c <^ 1 
when the star radius R is much smaUer than the light cylinder. As a result, one writes down the flrst six boundary conditions 
as 

e^(i?s,e) = 0, (31) 
Cj(Ps,e) = 0, (32) 
7^(Ps,e) = 7in, (33) 

i.e. £,^{Rs,9) — l/(27f„). According to all theories of particle generation near the neutron star surface (Ruderman Sutherland 
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1975, Arons Scharlemann 1979), ■jin < 10^ for secondary plasma. For this reason in what follows we consider in more details 
the case 

7in « ^, (34) 

when the additional acceleration of particles inside the fast magnetosonic surface takes place (see e.g. Beskin Kuznetsova 
Rafikov 1998). It is this case that can be realized for fast pulsars. Moreover, it has more general interest because the relation 
( p^ may be true also for AGNs. As to the case 7f„ 2> u corresponding to ordinary pulsars, the particle energy remains 
constant (7 — "fin) at any way up to the fast magnetosonic surface (see Bogovalov 1997 for details). 
Further, one can put 

S{R,,e) = 0, (35) 
ef{R,,6) = 0, (36) 
ri+{R,,9)-ri''(R,,9) = 0. (37) 

These conditions result from the relation cEs + QR^e^p x Bs = corresponding rigid rotation and perfect conductivity of the 
surface of a star. Finally, as will be shown in Sect. 3.2, the derivative dS/dr actually determines the phase of plasma oscillations 
only and plays no role in the global structure. Finally, the determination of the electric current and, say, the derivative df/dr 
depend on the problem under consideration. Indeed, as is well-known, the cold one-fluid MHD outflow contains two singular 
surfaces, Alfvenic and fast magnetosonic ones. It means that for the transonic flow two latter functions are to be determined 
from critical conditions (Heyvaerts 1996). In particular, the longitudinal electric current within this approach is not a free 
parameter. On the other hand, if the electric current is restricted by some physical reason, the flow cannot pass smoothly 
through the fast magnetosonic surface. In this case, which can be realized in the magnetosphere of radio pulsars (Beskin et al 
1983, Beskin & Malyshkin 1998), near the light surface \E\ = \B\ an effective particle acceleration may take place. Such an 
acceleration will be considered in Sect. 4. 



3 THE ELECTRON-POSITRON OUTFLOW 
3.1 The MHD Limit 

In the general case Eqns. ( |l9| ) - ( p^ ) have several integrals. Firstly, Eqns.(pl|), (pj]), and ( p6| ) result in 

2 (A -1/2 cos 0)7+ + (A +1/2 cos e)7- _ 1 , m 

^ tane 2crAsin6l a sin 61 sin 61 ' ^ ' 

where l{6) describe the disturbance of the electric current at the star surface by the equation I{R,6) = I a [sin"^ 6 + /(S)] . 
Expression ( p8| ) corresponds to conservation of the total energy flux along a magnetic field line. Furthermore, Eqns. (|2^) - 
( pij ) together with the boundary conditions (|35|), ( |3^ result in 

1+/ Ursine 
4Aa^ V c ^'^J^iXa 

' = ^^+4A^^ (l-^^O^iA^^- 

They correspond to conservation of the z-component of the angular momentum for both types of particles. It is necessary to 
stress that the complete nonlinearized system of equations contains no such simple invariants. 

As (tA ^ 1, we can neglect in Eqns. (^^-(^8|) their left-hand sides. In this approximation we have = ^~ i.e. 7" = 
7+ = 7, so that 

1 95 C sin 61 ^ c ^ „ , , _ 

' ^ '^ + 7^e«^0, (42) 



3 - e/-— 7in; (39) 



0,rsmO dr Qr^ 
and 

/ r2rsin6'^ \ , , 

7(1 —^^) = (43) 

Hence, within this approximation 

S = ef, (44) 
< - t^^^ + S-^(^-^-)- 

© 1999 RAS, MNRAS 000, [|-?? 



(Jn the particle acceleration near the light surface oj radio pulsars 5 

Substituting these expressions into and using Eqns. we obtain the following equation describing the disturbance 

of the magnetic surfaces 

e{l-x^ sin" e)^+e{l~x^ sin"" 9) ^ 4 i ^ ^1 - 2ea: sin' - 2e sin e cos + 2e(3 cos" 61-1)/ (46) 

+ ii^5^('"'^' - (7 - 7nO - - 2Asin' - ir) + ^ sine(C - ^,-) = 0, 

where a; = Qr/c. One can see that it actually coincides with the one-fluid MHD Eqns. (32), (52) from Beskin et al (1998), but 
contains the two last additional nonhydrodynamical terms. Nevertheless, as will be shown in the next subsection, at small 
distances r ^ rp where rp is the radius of the fast magnetosonic surface we have 

- A sin' e{^+ + ^ sine(C+ - C^) ~ 0, (47) 

so actually there is perfect agreement with the MHD approximation 

£(l-x'sin'e)f|+e(l-a;'sin'e)^^ (^^) - 2ex sin' - 2e sin cos (48) 
OX'' x^ oO \ sm 00 / ax aO 

As was shown earlier (Beskin et al 1998), to pass through the fast magnetosonic surface it's necessary to have 

\l\ < a-*''^ (49) 

Hence, within the fast magnetosonic surface r <C rp one can neglect terms containing S = ef and C^. Then, relations (^l|) and 
(^ result in 

7(1 -T sin e^^) = 7in, (50) 

6- = (51) 
a; sm 6' 

= 0. (52) 
Finally, using the definition 

we obtain for a 3> 7in for r <C rp 

7' = 7in + x' sin' 9, (54) 



v/7i'n + a^' sin' 8 - 7i„ 1 

= , ~ : — 7;, (55) 

ssin6i^72^ + a;2 sin' 6* 2;sm6' 

f, = \/7in + a:' sin' e ~ 7in ^ 1 .gg. 

i2 sin' 61 y'72„ + x^ sin' ^ ~ a;^ sin' ' 

in full agreement with the MHD results. 

Next, to determine the position of the fast magnetosonic surface rp, one can analyze the algebraic equations ( ^ ) and 
( ^l[ ) which give 

-^ + T^S ^7 - sin + = 0. (57) 

aO tan 6^ cr sm 6^ x 

Using now expressions ( ^3| ) and (|5^), one can find 

7' + asin'6» = 0, (58) 



27^ - 2a 



where 

Jf(r,6i) = 2cos6l5-sin6'|^. (59) 

Ou 

Equation (^8|) allows us to determine the position of the fast magnetosonic surface and the energy of particles. Indeed, 
determining the derivative rd^/dr, one can obtain 

dj ^ 1^ {rdK/dr - X-') 

dr 37 - a (2^- + a;-2) ' ^ ^ 

© 1999 RAS, MNRAS 000, [|-?? 



b V.b. Beskm and K.K. Kafikov 



As the fast magnetosonic surface is the X-point, both the nominator and denominator are to be equal to zero here. As a 
result, evaluating rdK/dr as K, we obtain 

5 ~ (j-^/^ (61) 
rp ~ (j'/^/Jl; (62) 
7(rF) = (ji/^ sin^/^ (63) 

where the last expression is exact. These equations coincide with those obtained within the MHD consideration. It is the 
self-consistent analysis when 5 — ef, and hence K depends on the radius r that resuhs in the finite value for the fast 
magnetosonic radius rp. On the other hand, in a given monopole magnetic field, when ef does not depend on the radius, the 
critical conditions result in rp ^ oo for a cold outflow (Michel, 1969, Li et al 1992). 
Near the fast magnetosonic surface r ~ o^^^Ri^ the MHD solution gives 

7 ~ (64) 
ef ~ o--'/\ (65) 



Hence, Eqns. (|53|), (pq), and ( |56D result in 

~ (66) 

Cfl ~ (67) 
~ (68) 

As we see, the ^-component of the velocity plays no role in the determination of the 7. 

However, analyzing the left-hand sides of the Eqns. (p^)-(p^) one can evaluate the additional (nonhydrodynamic) varia- 
tions of the velocity components 

AC^ ~ A^V-'*''^ (69) 
AC± ~ A-V-^/^ (70) 
ACj ~ A-V-\ (71) 

Hence, for nonhydrodynamic velocities A^,'^^ <C and A^J <C to be small, it is necessary to have a large magnetization 
parameter a ^ 1 only. On the other hand, A^^/^g ~ A~^. In other words, for a highly magnetized plasma a 2> 1 even outside 
the fast magnetosonic surface the velocity components (and, hence, the particle energy) can be considered hydrodynamically. 
The difference ~ A~^ appears in the 9 component only, but it does not affect the particle energy. Finally, one can obtain from 
(H), (0) that 

~ A-^a-^/^ (72) 

To put it differently, at large distances the nonhydrodynamical terms are much smaller than hydrodynamical ones. 

As a result, at large distances where, according to (p9[)-(|4C|), one can neglect the toroidal component we obtain 

S = £/, (73) 

C = -^5-cr~'^7- (74) 
^ tan e sine ' ^ ' 



On the other hand, Eqn. ( |23[ ) gives 

C=|^+sineC.. (75) 
Together with (ji^) one can obtain for r ^ rp 

7 = cr(^2cose»e/-esin6i|^^ , (76) 

which coincides with the MHD solution. Finally, using Eqns. ([l9|), (|2^), and neglecting the nonhydrodynamic term 4A(^^ — ^^), 
one can flnd 

4r {^'%) -'^'^'^^ - ^-^1^^ + = °- 

Together with ( [t^) this equation in the limit r ^ rp coincides with the asymptotic version of the trans-field equation 
(Tomimatsu 1994, Beskin et al 1998) 

where g{e) = K{e)/sm'^ 6, and 
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D + i= , . 4. g"'W« 1- (79) 

sm 

In this limit, none of the terms containing and plays role in the asymptotic trans-field equation. Ifence, it is not 
necessary to consider the effect of the nonhydrodynamical term 4A(^^ — either. 



3.2 Plasma Oscillations 

In the intermediate region r <C rp Eqn. ( fF^ ) cannot be used. The point is that in the limit A 3> 1 the important role in Eqns. 
( p^ and (|2^) is played by the nonhydrodynamic terms (^) corresponding to different velocities of two components. As a 
result, the full version of Eqn. (0) has the form 

I- f^'^) -4coseC,-sin0§ + ^-^(C^sine) + 2A(C+-Cr) =0. (80) 
or \ or J OO xsmO oO 



Indeed, one can see from equations ( |l9| ) and ( |20[ ) that near the origin a; = i?s in the case 7j"[^ = 7;~ (and for the small variation 
of the current ( ~ cr"*''^ which is necessary, as was already stressed, to pass through a fast magnetosonic surface) the density 
variation on the surface is large enough: (r;^ — 77") ~ 7j~^ 3> C- Hence, the derivative d^S/dr^ here is of the order of 7;"^. On 
the other hand, according to (p^), the derivative ed^f/dr^ is x'^ times smaller. This means that in the two-component system 
the longitudinal electric field is to appear resulting in a redistribution of the particle energy. Clearly, such a redistribution 
is impossible for the charge-separated outflow. In other words, for a finite particle energy a one-component plasma cannot 
maintain simultaneously both the Goldreich charge and Goldreich current density (^) . In a two-component system with A ^ 1 
it can be realized by a small redistribution of particle energy (Ruderman & Sutherland 1975, Arons & Scharlemann 1989). 

For simplicity, let us consider only small distances a:: ^ 1. In this case one can neglect the changes of the magnetic 
surfaces. Using now (M) and (|2(l), we have 



7+ = 7in-4Acr<5; (81) 
7" = 7in+4Acr<5. (82) 
Finally, taking into account that ^9 and are small here, one can obtain from ( po| ) 



2d^s . „ as , 1 d f . „dS 



+ 2r-— + 



(sin.|^)+A.^^, (83) 



Q^2 Qj. gjjj 6 89 \ 89 J 7; 
where 

A = 16^ » 1. (84) 
Eqn. (^3|) has a solution 

5 = So + r^^^'^ [Ci sm{fi\nr) + C2 cos(/ilnr-)] cos 5, (85) 

where 

7in cos 9 
^»"T6A5^' 

and ^1 ~ ^/A. As we see, Eqn. (^) describes plasma oscillations similar to those considered by Shibata (1997) for charge- 
separated flow. The decrease of oscillations results from a more accurate consideration of the Laplace operator in a 3D 
space. 

One can easily check that the additional potential So is small, and it is not necessary to add it in (jsij) and (^9|)-(^o|). 
Moreover, the nonhydrodynamic disturbance A^^ (as well as A7) is also small, A^r/Cr ~ A~^. Hence, as was already stressed, 
the boundary condition 05 /Or (determining together with (^) the coefficients Ci and C2) does not affect the general structure 
of the flow. On the other hand, the presence of an additional electric potential So results in a full compensation of the last 
term in dTgl) 



2A(?+ ~(.r)- cos 9^r ^ 0. (87) 
Next, as ef ^ a^^^'' for r <^ rp, a similar expression can be written for the (^-components as well 

2A(?J-C^)-cos6»e^RiO. (88) 

Expressions (|^ - ( p8[ ) must hold for the whole region r < rp. In this case, the final version of Eqn. (|o|) in the internal region 
r <C rp can be rewritten as 

l-fr^l^) -2coseC,-sine^C, + ^-4(C^sin0) = O. (89) 
8r \ 8r J 09 xsm 9 89 
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As 5 ~ e/ <C (7 '^^'^ for r <C rp, and ~ 7q ^ ^ 5, the first term in (^) can be omitted. As a result, the solution of 
Eqn. (H) coincides exactly with the MHD expression, i.e. 7^ = ■yf^ + x'^ sin^ 6 Finally, using (|^), (Q, and (||)-(|6|), 

one can easily check that the nonhydrodynamical terms in the trans-field equation ( |48| ) do actually vanish. 



4 THE BOUNDARY LAYER 

Let us now consider the case when the longitudinal electric current I(R, 9) in the magnetosphere of radio pulsars is too small 
(i.e. the disturbance l{9) is too large) for the flow to pass smoothly through the fast magnetosonic surface. First of all, it can 
be realized when the electric current is much smaller than the Goldreich one. This possibility was already discussed within 
the Ruderman-Sutherland model of the internal gap (Beskin et al 1983, Beskin & Malyshkin 1998). But it may take place 
in the Arons model (Arons & Scharlemann 1979) as well. Indeed, within this model the electric current is determined by 
the gap structure. Hence, in general case this current does not correspond to the critical condition at the fast magnetosonic 
surface. In particular, it may be smaller than the critical current (of course, the separate consideration is necessary to check 
this statement). 

For simplicity let us consider the case l{9) = hsiv? 9. Neglecting now the last terms oc a^^ in the trans-field equation 
(PI), we obtain 

,(1_, sm 9)—+eil^a: sm g) — - -j 

-2ex sin^ - 2e sin 6* cos 6*^ + 2e(3 cos^ 61 - 1)/ + 4/i sin^ 9 cos 9^ 0, (90) 
ox o9 

which actually coincides with the force-free equation (Beskin et al 1998). This equation has an exact analytical solution 

e f = hx^ sirl^ 9 cos 9 . (91) 

For h < (when the electric current is smaller than the Goldreich one) this solution results in the appearance of the light 
surface \E\ = \B\ at the finite distance 

As we see, for l{9) = hsin^ 9 this surface has the form of a cylinder. It is important that the disturbance of magnetic surfaces 
ef ~ (I'll)^^^ remains small here. 

Comparing now the position of the light surface (^2|) with that of the fast magnetosonic surface (|6^, one can find that 
the light surface is located inside the fast magnetosonic one if 

a-"^^ « \h\ « 1, (93) 

which is opposite to (^). One can check that the condition (^^ just allows to neglect the non force-free term in Eqn. (^). 
Using now the solution (|9l[) and the MHD condition 5 = ef, one can find from (BSl) 



27^ - 2(7 (^hx'^ sin" 9 + 7^ + 0- sin^ 9 = 0. (94) 

This equation shows that near the force-free boundary Xff = (2|ft|)^^/* ( [j^ ) 

7 = sin^/3 9 - ^M^!^V3 ^-^^/s e^xn~xsm9, (95) 
v3 

where 

if Q 1 

(96) 



xq = 



(2|/i|)V4 



4(2!/ii)i/2 (crsin^ 61)2/3 



(see Fig. |l]). Hence, the real solution is absent for a; sin 6' > xo. Here 7(2:0) = ct^''"' sin^''"' 6', the condition a"*^^ <^ \h\ resulting 
in zoc < rp, and the last term in (^) being small. 

Although the energy of particles at the limiting point is finite, the derivative d'y /dm moves to infiruty. Hence, near the 
light surface the left-hand sides in the Eqns. ( p^ - (^sj) are to be taken into consideration. Since in our case the light surface 
has the form of a cylinder, one can move to derivatives perpendicular to the boundary layer only by 

d/dr ^ sin 9d/dvj; (97) 
d/d9 ^ zocos9d/dvj. (98) 

As a result, ^ can be eliminated from (^) and (|2l|). Together with (^) they give the equation for 5 (see (104)). Next, the 
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hotindaiy layer 
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Figure 1. The behavior of the Lorentz factor in the case tr <C <C 1. One can see that the one— fluid MHD solution @ exists 
for 7 < (T^/^ only. But in the two— fluid approximation in the narrow layer Aro = tuc/A the particle energy increases up to the value ~ a 
corresponding to the full conversion of the electromagnetic energy to the energy of particles. 



invariants (p3 ) and (Uol) can be used to define ^J: 



X sin 6 
1 

X sin S 



1 + 



4Aa(<5 - £/) 



1 



7+ 

4Acr((5 - £/) 



7" 



Furthermore, one can define 

As to the energy integral (^sj), it determines the variation of the current (. Now it can be rewritten as 



2 _g (7+ + 7~) 



tan 9 2(7 sin 9 
Finally, Eqns. (Ili) - (Ei) look like 



tj?;,-; — TT = 2 sin ^ cos ^ 



-2sin''e' 

-2sin^6' 

„+ 



(99) 
(100) 

(101) 
(102) 

(103) 



A-icose) (A+icos^les 



7 



7 



(Ce-7" 



-4Aa 



(^A - i cos ^) - (-^ + ^ cos 5; 

(a cose) (A+icose) 

cos 6* d(5 . sine^+ 

^5 — ^-5- - ^=^-^3 sme^,. H 

icr sm 9 sm dro a;o 

7^ + 7~ cos 6 d5 . sinS,_ 

- — ^ - ^c^-^-: sme^^ + ^ 

2(7 smt? smt? dm a;o 

d + / dS . 4. 

roc-^7 = 4A(7 -OTc- sm 9t,g 

azu V dec 



(104) 
(105) 
(106) 

(107) 
(108) 
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d _ 

avj 



-4Ao- i^- 



dS 

' dzu 



sin O^^g 



where we neglected the terms (x S/r in (106) and (107). 

Comparing the leading terms, we have inside the layer A-w/Ri^ ~ A^^ 



1,1/2 

h^' a, 

where he = \h\. Then the leading terms in 
1 1 



7 

AS 



(103) for Aro > A"^i?L are 



X sin 6 
1 



Xo 

_ 1 

X sin 9 Xo ' 



2?+ = (C)' + (a)', 

2^.- = (5;f + (C9")', 

(7+ +7") 



c = - 



2a sin ( 



(109) 



(110) 
(111) 
(112) 
(113) 

(114) 
(115) 

(116) 
(117) 

(118) 



where xq = zuc /Rl = (2|/i|)-i/*. Hence, one can totally neglect ef and 5 in ( |99| ) - ( |lO(]| ), so Eqns. ( |l04| ) - ( |l09[ ) in the region 
Auj > A~^iiL can be rewritten as 



^2 d'5 



= 4Acr 



' dvj 



(Ce"7')=-4Aa - 



7^+7 
2(7 sin 6 

+ 7" 



7 



2(T sin ^ 



cos dS 
' sin 61 dzu 

cos 6 dS 
' sin 6 dzu 



dzu ^ ' 
d 

azu 





(119) 


- sin6'5+^ , 


(120) 


— sin 6^r^ , 


(121) 


4Xa sin 6(g , 


(122) 


4Acr sin 9^g , 


(123) 



with all the terms in the right-hand sides of (|l20|) an d (121) being of the same order of magnitude. 

As a result, the nonlinear equations (11!:) - (123) and (105) give the following simple asymptotic solution 

7* = 4smea{Xlf, 

=F2sin6'A;, 
4 . 



AS 

A(e/) 
C 



sin 6 cos 6X 



cos6lA"^(A^)^ 
'(AO', 



~4sin6'(A0 



(124) 
(125) 

(126) 

(127) 
(128) 



where now I = Azu/zuc- It is important that the last expressions are correct for arbitrary relations between 7?-; and a. As 
we can see, in the narrow layer Azu — zu^/X the particle energy increases up to the value ~ a which corresponds to the full 
conversion of the electromagnetic energy to the energy of particles. For this reason we have here \(\ ~ 1, which just means the 
diminishing of the toroidal magnetic field determining the flux of the electromagnetic energy. On the other hand, the variation 
of the electric potential remains small S ~ A~^, to say nothing about the variation of the magnetic surfaces Aef ~ A"'^. These 
results coincide exactly with our previous evaluations (Beskin et al 1983) allowing us to neglect variations of the electric 
potential and the poloidal magnetic structure in the ID cylindrical case. 

The latter result has a simple physical explanation. Indeed, the diminishing of the toroidal magnetic field is connected 
with the ^-component of the electric current which is produced by all the particles moving in opposite directions. On the 
other hand, the change of the electric potential depends on the small difference between the electron and positron densities. 
As a result, according to (127) and (128), the change of the toroidal magnetic field is just A times larger than the change 
of the electric potential. Unfortunately, it is impossible to consider this region more thoroughly because for A? ~ 1 we have 
5^ ~ 1 and ~ 1, i.e. the linear approximation (^) - (j28|) itself becomes incorrect. 

It is necessary to stress as well that we do not include into consideration the radiation reaction force 
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(rad) 



3m2c4 



[{Ey~B.f + {E-.^Byf] 



(129) 



which can be important for large enough particle energy. Comparing (129) with appropriate terms in (120) - (125) one can 
conclude that the radiation force can be neglected for o < (Jcr, where 

i lo^ (130) 



and re — jnu? 



<3 X 10"'' 



- classical electron radius. This relation can be rewritten in the form 

-3/7-, 2/7 



which gives 



p > o.oeSi^ ' A 



3/7 > -2/7 



(131) 



(132) 



Hence, for most radio pulsars the radiation force indeed can be neglected. As to the pulsars with a > (Jcr, it is clear that for 
7 > CTcr the radiation force becomes larger than the electromagnetic one and strongly inhibits any further acceleration. As a 
result, we can evaluate the maximum gamma-factor which can be reached during the acceleration as 



7ii 



10° 



(133) 



5 DISCUSSION 

Thus, on a simple example it was demonstrated that for real physical parameters of the magnetosphere of radio pulsars 
(a 3> 1 and A S> 1) the one-fluid MHD approximation remains true in the whole region within the light surface \E\ — \B\. 
On the other hand, it was shown that in a more realistic 2D case the main properties of the boundary layer near the light 
surface existing for small enough longitudinal currents I < Iqj (effective energy transformation from electromagnetic field to 
particles, current closure in this region, smallness of the disturbance of electric potential and poloidal magnetic field) remain 
the same as in the ID case considered previously (Beskin et al 1983). 

It is necessary to stress the main astrophysical consequences of our results. First of all, the presence of such a boundary 
layer explains the effective energy transformation of electromagnetic energy into the energy of particles. As was already 
stressed, now the existence of such an acceleration is confirmed by observations of close binaries containing radio pulsars (as 
to the particle acceleration far from a neutron star, see e.g. Kennel & Coroniti 1984, Hoshino et al 1992, Gallant & Arons 
1994). Simultaneously, it allows us to understand the current closure in the pulsar magnetosphere. Finally, particle acceleration 
results in the additional mechanism of high-energy radiation from the boundary of the magnetosphere (for more details see 
Beskin et al 1993). 

Nevertheless, it is clear that the results obtained do not solve the whole pulsar wind problem. Indeed, as in the cylindrical 
case, it is impossible to describe the particle motion outside the light surface. The point is that, as one can see directly from 



Eqn. (126), for a complete conversion of electromagnetic energy into the energy of particles it is enough for them to pass 
only A~^ of the total potential drop between pulsar magnetosphere and infinity. It means that the electron-positron wind 
propagating to infinity has to pass the potential drop which is much larger than their energy. It is possible only in the presence 
of electromagnetic waves even in an axisymmetric magnetosphere which is stationary near the origin. Clearly, such a flow 
cannot be considered even within the two-fluid approximation. In our opinion, it is only a numerical consideration that can 
solve the problem completely and determine, in particular, the energy spectrum of particles and the structure of the pulsar 
wind. Unfortunately, up to now such numerical calculations are absent. 
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